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ABSTRACT 


The  heat  transfer  characteristics  and  performance  capabilities  of  a  liquid  flow¬ 
through  module  for  electronics  cooling  were  investigated  experimentally  using  the 
Standard  Electronics  Module,  format  "E"  (SEM-E)  size  module.  Brayco  Micronic 
889,  a  dielectric  polyalphaolefin,  was  tested  as  the  coolant.  One  surface  of  the 
flow-through  module  was  populated  with  six  etched  foil  heaters  placed  on  the 
module  over  the  fluid  flow  path.  With  thermocouples  affixed  to  the  surface  and 
placed  in  the  inlet  and  outlet  tubing,  data  was  gathered  and  analyzed  to  quantify 
module  effectiveness  in  dissipating  heat  over  a  variety  of  power  input  settings  and 
coolant  flow  rates. 
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I. 


INTRODUCTION 


A.  ELECTRONICS  COOLING 

Modern  advances  in  microelectronic  technology  have 
steadily  increased  the  heat  dissipation  levels  of  integrated 
circuit  chips  and  mandate  improved  thermal  management 
techniques .  Heat  fluxes  at  the  surface  of  a  chip  may  reach 
200  W/cm^  and  as  efforts  to  miniaturize  electronics  continue 
even  greater  fluxes  are  expected.  Equipment  reliability 
becomes  a  major  concern  as  components  are  subjected  to  higher 
temperatures.  One  study  showed  that  temperature  fluctuations 
greater  than  20°C  about  an  average  temperature  resulted  in 
increasing  the  failure  rate  eightfold  (Hilbert  and  Kube, 
1969)  .  Failure  may  result  from  thermal  fracture  of  the 
usually  ceramic  or  metallic  substrate,  plastic  flow  of  printed 
circuits,  etc.,  due  to  elevated  temperatures,  even  without  the 
presence  of  such  fluctuations.  In  order  to  hold  failure  rates 
in  the  vicinity  of  0.5  to  2.0%  per  thousand  hours  of 
operation,  component  junction  temperatures  of  over  100°C  have 
been  commonly  accepted  in  the  past,  yet  by  further  decreasing 
temperatures,  one  can  increase  expected  life  and  decrease 
maintenance  costs  as  well  as  associated  down  time.  (Kraus  and 
Bar-Cohen,  1983) . 

Natural  convection  and  radiation  cooling  relying  upon  the 
ambient  air  is  insufficient  to  deal  with  the  rapidly 
increasing  fluxes.  Space  limitations,  which  often  discourage 
the  use  of  conventional  forced  convection  air  cooling,  provide 
an  impetus  to  the  development  of  alternative  methods  of 
cooling.  Better  thermal  performance  can  also  pave  the  way  for 
weight  savings  as  the  density  of  circuitry  might  be  allowed  to 
increase  without  increasing  component  housing  dimensions. 
Such  improvements  have  readily  visible  advantages  for  ship  and 
aircraft  design  where  space  is  at  a  premium. 

Much  research  is  currently  underway  in  the  field  of 
electronics  cooling  as  attempts  are  made  to  address  these  many 
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concerns,  with  the  need  to  bring  a  capable  cooling  medium 
closer  to  the  critical  components  at  its  focus.  This  study- 
focuses  on  the  use  of  liquid  flow-through  modules  and  their 
effectiveness  in  heat  removal  using  a  polyalphaolef in  coolant. 

1.  Liquid  Cooling 

Natural  convection  cooling  in  air  is  effective  for  only 
relatively  small  heat  fluxes,  yet  forced  convection  air 
cooling  can  only  increase  heat  dissipation  to  the  vicinity  of 
1  W/cm^  if  allowing  a  100°C  temperature  difference  (Kraus  and 
Bar- Cohen,  1983) .  Liquid  cooling  has  thus  emerged  as  an 
attractive  alternative  method  of  cooling  with  direct  immersion 
and  conduction  cooling  at  the  center  of  much  research. 
Because  of  the  greater  heat  transfer  coefficients  associated 
with  liquids,  their  use  in  cooling  allows  for  lengthier 
emergency  operation  of  vital  equipment  in  the  event  of 
interrupted  flow.  Some  disadvantages  do  exist  with  liquid 
cooling,  such  as  the  added  system  weight  and  additional  costs 
associated  with  the  fluid  flow  system.  There  may  also  be  the 
potential  for  damage  to  electronic  components  due  to  system 
leaks.  Lightweight  dielectric  coolants,  such  as  the  Brayco 
Micronic  889  used  in  this  study,  are  therefore  promising 
solutions  to  the  thermal  challenges  confronted  in  the  design 
of  electronics  cooling  systems.  Brayco  Micronic  889  is  a 
polyalphaolef  in  made  for  use  in  closed  loop  systems  for  ground 
based  and  airborne  electronics  cooling.  It  lays  claim  to  high 
specific  heats  and  to  compatibility  with  typical  construction 
metals.  (Castrol,  1994). 

2 .  Finned  Heat  Exchangers 

Extended  surfaces,  or  fins,  have  become  a  commonplace  aid 
in  convective  heat  transfer.  Though  abundant  research  has 
sought  to  improve  convective  heat  transfer  effects  using  fins, 
most  past  studies  have  dealt  with  the  use  of  air  and  other 
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gases.  The  usefulness  of  such  data  is  very  limited  when 
evaluating  the  performance  of  liquid  coolants.  Because  of  the 
higher  Prandtl  numbers  and  thermal  conductivities  generally 
associated  with  liquid  coolants,  the  methods  of  enhancing  heat 
transfer  are  slightly  different  than  those  methods  used  with 
gases.  The  differences  between  liquids  and  gases  in  the 
behavior  of  properties  such  as  density  and  specific  heat  also 
necessitate  additional  research.  With  liquids,  rather  than 
simply  increasing  the  "wetted"  surface  area  through  branching 
mazes  of  fins  and  narrowing  the  resultant  channels,  the 
convection  coefficient  may  best  be  elevated  by  the 
interruption  of  the  developing  thermal  boundary  layers.  Heat 
exchangers  such  as  the  flow-through  module  use  offset  plate 
fins,  similar  to  those  shown  in  Figure  1,  to  enhance  heat 
transfer  through  the  periodic  starting  and  development  of 
laminar  thermal  boundary  layers  over  flow  channels  formed  by 
the  fins  and  their  dissipation  in  the  fin  wakes.  This  flow  is 
also  coupled  with  an  increase  in  the  pressure  drop  due  to  the 
friction  and  form  drag  associated  with  the  fins. 

3 .  Flow-Through  Modules 

Offset  plate  fins  are  the  core  of  the  design  of  a  liquid 
Flow-Through  Module  (FTM)  such  as  the  Standard  Electronics 
Module,  Format  E  (SEM-E)  size  module.  Flow-through  modules 
are  manufactured  through  vacuum  brazing,  which  allows  for  the 
creation  of  very  narrow  channels  and  fins  in  a  thin  module. 
The  module  tested  is  0.100  inches  thick,  with  an  0.02  inch 
Aluminum  6061  faceplate.  The  fins  within  the  module  are  made 
from  an  Aluminum  3003  alloy  and  are  0.06  inches  in  height, 
0.125  inches  long,  and  a  mere  0.006  inches  thick.  The  fin 
density  is  twenty  fins  per  inch.  (Buechler  and  Brough,  1993)  . 

Electronic  components  placed  on  the  surface  of  a  module 
are  cooled  by  the  coolant  flowing  within  the  module.  Both 
surfaces  of  the  module  could  be  used  in  this  manner  with  the 
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added  benefit  of  the  structural  support  to  the  circuitry- 
provided  by  the  plate .  The  thin  module  structure  is  very 
suitable  for  stacking.  Several  modules  can  be  housed  in  a 
small  enclosure  with  rails  to  guide  insertion  and  maintain 
some  separation.  The  module  manufacturing  process  allows  for 
the  internal  directing  of  coolant  to  areas  of  particular 
thermal  concern  on  the  surface  through  the  flexibility  of  flow 
path  design.  The  FTM  tested,  shown  in  Figure  2,  possesses 
three  main  flow  paths  filled  with  fins  which  channel  the 
coolant  along  a  winding  path  from  entrance  to  exit  (Figure  3)  . 
The  flow-through  module  is  currently  under  evaluation  by  the 
military  for  increased  application. 


B .  OBJECTIVES 

This  study  examined  the  effectiveness  of  the  SEM-E  sized 
flow-through  module  using  a  synthetic  dielectric 
polyalphaolef in  coolant  known  as  Brayco  Micronic  889, 
manufactured  by  Castrol.  Several  different  flow  rates  were 
utilized  while  varying  the  heat  flux  supplied  to  the  surface 
by  etched  foil  heaters .  Temperature  data  for  the  module 
surface  was  collected  in  order  to  quantify  the  effectiveness 
of  the  module  and  provide  the  groundwork  for  potential  further 
studies  into  the  analytical  modelling  of  the  convection 
process  occurring  within  the  module.  Data  on  inlet  and  outlet 
temperatures  was  also  collected  for  analysis  of  the  heat 
dissipation  capabilities  of  the  module  over  a  variety  of  power 
inputs  and  flow  levels. 
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Figure  1.  Offset  plate  fins  (from  Herold,  Srindar, 

and  Hu,  1992)  . 
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Module  Tiiickiiess:  0.1  incli 
Faceplate  Tkicknessr  0.02  incli 
Faceplate  Material:  6061  Alum. 

Figure  3.  Fluid  flow  path  within  FTM  (from  Buechler  and 

Brough,  1993)  . 
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II.  EXPERIMENTAL  APPARATUS 


A.  SUPPORT  SYSTEM  ASSEMBLY 

The  experimental  apparatus  consisted  of  the  flow-through 
module  assembly  and  the  support  system  as  drawn  in  Figure  4 . 
A  description  of  each  sub-system  and  its  components  follows. 

1.  Fluid  Circulation  System 

An  Endocal  RTE-5  refrigerated  circulating  bath  provided 
the  heat  sink  for  the  coolant,  reducing  its  temperature  to 
approximately  ambient  room  temperature  prior  to  entering  the 
pump.  The  pump  was  a  Masterflex  variable  speed  drive  with  an 
attached  pump  head,  providing  a  flow  of  up  to  one  liter  per 
minute  through  flexible  Tygon  tubing.  A  Fischer  and  Porter 
flow  meter  was  used  to  determine  the  flow  rate.  Calibration 
was  accomplished  by  timing  the  collection  of  a  volume  of 
coolant  over  several  meter  readings .  The  maximum  uncertainty 
in  the  measured  flow  in  the  range  of  testing  was  2.8%. 
Calibration  data  was  then  fitted  to  a  fifth  order  polynomial 
using  TableCurve  software.  The  polynomial  was  used  to 
determine  flowmeter  settings  for  a  desired  rate  of  flow  during 
experiments.  Appendix  B  includes  flow  meter  uncertainty 
calculations  and  the  polynomial  curve  fit. 

Brass  tubing  was  threaded  to  fit  the  inlet  and  exit  ports 
of  the  module  in  order  to  provide  an  attachment  point  for 
tubing.  The  brass  tubes  were  inserted  into  the  Tygon  tubing 
and  clamps  were  slipped  over  the  connection  and  tightened. 
Inlet  and  exit  tubing  was  fed  into  Teflon  T-fittings,  each  of 
which  was  modified  to  accommodate  a  thermocouple.  Upon 
exiting  the  module  and  passing  through  the  T-fitting,  the 
coolant  returned  to  the  refrigerated  bath.  Interconnections 
between  components  were  made  with  plastic  and  Tygon  tubing  and 
sealed  with  an  acrylic  adhesive.  Clamps  were  affixed  to 
connections  to  maintain  system  integrity  at  high  pressures. 


9 


2.  Power  Distribution  System 

Power  to  the  heaters  was  supplied  by  a  Hewlett  Packard 
0-20  V,  0-10  A  DC  power  supply.  The  six  heaters,  wired  in 
parallel,  were  7.6  mm  (0.3  inch)  by  38.1  mm  (1.5  inches)  Minco 
Kapton  etched  foil  heaters.  Two  heaters  were  applied  with 
thermally  conductive  epoxy  adhesive  (Omegabond  101)  along  each 
of  the  three  module  flow  channels  at  a  spacing  of  5  cm 
(2  inches) .  Flow  channels  were  roughly  determined  by  flow 
visualization  utilizing  a  transparent  polymer  sheet  of 
Hallcrest  thermochromic  liquid  crystals.  Each  row  of  heaters 
was  separated  by  1.27  cm  (0.5  inch).  Each  heater  had  an 
effective  heated  area  of  1.94  cm^  (0.3  in^)  and  a  resistance 
of  20.7  ohms  +  10%.  Figure  5  shows  the  actual  module  once 
heaters  were  attached.  The  side  of  the  module  to  which 
heaters  were  affixed  was  arbitrarily  designated  as  side  "A." 

Thermally  stable  precision  resistors  with  resistances  of 
0.487  ohms  +  1%  were  connected  in  series  with  each  heater. 
Each  resistor  was  tested  using  a  Rosemount  920A  Commutating 
Bridge  to  establish  its  resistance  (R)  to  a  fourth  decimal 
place. 

An  HP  3852A  Data  Acquisition/Control  Unit  with  an  HP 
44701A  digital  voltmeter  was  used  to  determine  voltages.  The 
voltage  drop  across  a  precision  resistors  (Vj.)  was  first 
measured  utilizing  an  HP  44705A  multiplexer  and  then  current 
to  the  heater  was  determined  using: 
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Voltage  was  also  measured  across  each  heater  using  the  HP 
44705A  multiplexer  and  HP  44701A  voltmeter.  Power  supplied  to 
each  heater  was  calculated  as ; 

Ph  =  V,  Ih  (2) 

3 .  Data  Acquisition  System 

A  personal  computer  fitted  with  a  National  Instruments 
GPIB-PCIIA  IEEE-488  interface  board  was  used  to  control  an 
HP3852A  Data  Acquisition/Control  Unit  for  measurement  of 
voltages  and  temperatures.  Data  runs  were  performed  utilizing 
a  program  in  QBASIC.  An  HP  44708A  20  channel  relay 
multiplexer  for  thermocouple  compensation  provided  the 
interface  between  the  HP  44701A  voltmeter  and  the 
thermocouples  for  temperature  readings.  Fifteen  thermocouples 
were  on  the  module  itself.  In  addition  one  thermocouple  was 
in  the  inlet  and  one  was  in  the  outlet  tubing.  Two 
thermocouples  were  in  the  insulation  layers  above  the  module 
to  estimate  heat  loss.  An  additional  thermocouple  was  on  the 
surface  of  the  outer  layer  of  insulation,  thus  utilizing  all 
twenty  channels  of  the  multiplexer.  Figure  6  shows  a  layout 
of  the  thermocouple  locations  as  well  as  the  heater  locations. 
Thermocouples  adjacent  to  the  heaters  were  placed  as  close  to 
the  heaters  as  possible,  while  the  other  module  thermocouples 
were  placed  midway  between  the  heaters . 

Thermocouple  extensions  were  fed  into  channel  numbers  0 
through  19  of  the  HP  44708A  in  the  data  acquisition  system. 
Heater  voltages  were  measured  by  channels  0  through  5,  which 
corresponded  to  heaters  A  through  F  respectively.  Precision 
resistor  voltages  were  measured  utilizing  channels  6  through 
11  of  the  HP  44705A  relay  multiplexer. 
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B.  TEST  SECTION  ASSEMBLY 


1 .  Module  Housing 

A  plexiglass  enclosure  was  made  to  support  the  FTM  and 
insulation  in  a  horizontal  position,  thus  preventing  natural 
convection  effects  which  might  occur  in  an  upright  or  inclined 
system.  The  sides  of  the  structure  were  made  from  0.5  inch 
(1.25  cm)  thick  plexiglass  and  the  back  was  made  from  0.125 
inch  (0.32  cm)  thick  plexiglass.  The  module  and  insulation 
rested  upon  a  horizontal  0.25  inch  (0.64  cm)  plexiglass  piece 
which  was  secured  to  the  sides  by  screws.  Insulation  was  one 
layer  of  0.75  inch  (1.91  cm)  black  foam  insulation  below  and 
two  layers  above  the  module  with  a  layer  of  insulation  along 
the  sides  and  back  of  the  module  as  well.  The  enclosure  of 
the  module  was  open  in  the  front  to  allow  for  protruding 
wiring  and  tubing. 

2 .  Instrumentation 

Copper-constantan  (T-type)  thermocouples  were  affixed  to 
the  upper  surface  of  the  module  as  viewed  in  the  housing. 
Flat,  5  mil  thick,  foil  type,  cement -on  thermocouples  were 
used  to  allow  placement  very  close  to  the  heaters. 
Thermocouple  leads  were  separated  using  Kapton  tape.  Copper 
leads  were  soldered  to  thermocouple  extension  wire  and  the 
constantan  leads  were  resistance  welded  to  the  extension  wire. 
These  connections  were  also  wrapped  in  Kapton  tape  to 
guarantee  their  separation.  Thermally  conductive  paste  was 
applied  to  the  thermocouple  junction  at  the  application  point 
to  minimize  any  air  resistance  and  Kapton  tape  was  placed  over 
each  t he rmio couple  to  prevent  dislocation. 

Coolant  temperature  was  measured  using  T-type 
thermocouple  probes  in  the  inlet  and  outlet  tubing  from  the 
module.  The  probes  were  made  from  8  mil  thermocouple  wiring 
arc  welded  to  form  a  "bead"  junction.  The  probes  were  fed 
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through  a  small  hole  drilled  in  a  plastic  slug  and  inserted 
into  a  Teflon  T-fitting.  The  beads  were  just  below  the  center 
of  the  flow  path  in  order  to  achieve  a  more  "average"  value 
for  fluid  temperature.  Small  wire  screen  filters  were  made 
for  the  T- fittings  to  prevent  any  small  particles  from 
entering  the  flow.  One  filter  was  designed  to  rest  on  the 
inner  shoulder  of  the  T-fitting,  upstream  of  the  inlet 
thermocouple  probe.  A  filter  was  placed  upstream  of  the 
thermocouple  probe  in  the  outlet  T-fitting  in  a  similar 
fashion  in  order  to  promote  fluid  mixing  prior  to  the 
thermocouple  measurement . 

Thermocouples  placed  above  the  module  for  determining 
losses  through  the  insulation  were  made  from  24  gage 
copper- constant an  extension  wire  arc  welded  to  form  bead 
junctions.  These  thermocouples  were  placed  along  an  imaginary 
vertical  line  extending  above  thermocouple  13  as  shown  in 
Figure  6.  This  allowed  for  a  direct  flow  path  to  be  examined 
for  estimating  conduction  losses  through  the  insulation. 
Electrical  tape  was  used  to  hold  the  thermocouples  in  place  on 
the  insulation. 

All  thermocouples  were  initially  calibrated  in  an  ice 
bath  prepared  from  distilled  water.  Using  the  data  collection 
system,  no  thermocouple  reading  differed  by  more  than  0.11°C 
from  0 . 0°C  when  placed  in  the  ice  bath.  Because  of  the  desire 
for  the  highest  level  of  accuracy  in  the  determination  of 
fluxes,  the  inlet  and  outlet  thermocouples  were  additionally 
calibrated  using  a  Rosemount  920A  commutating  bridge  and  a 
constant  temperature  ethylene  glycol  bath.  Both  thermocouples 
were  placed  in  the  bath  and  tandem  readings  collected  over  a 
temperature  range  of  17°C  to  46°C.  These  values  were  compared 
with  the  reference  value  obtained  from  a  platinum  resistance 
thermometer  also  in  the  bath.  The  results  indicated  that  the 
the  two  thermocouples  produced  readings  within  0 . 03°C  of  each 
other,  although  they  varied  by  as  much  as  0.13°C  from  the 
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reference.  Table  26  in  Appendix  B  shows  the  readings 
collected  during  calibration.  These  figures  are  used  in 
determining  overall  uncertainty  as  discussed  in  Appendix  B. 
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Figure  4.  Diagram 


experimental  setup. 


INLET  OUTLET 

Heater  locations  are  indicated  by  lettered  rectangles. 

Spacing  is  approximately  2"  (5  cm)  horizontally  and  0.5"  (1 .25  cm)  vertically. 
Thermocouple  locations  are  indicated  by  numbers. 
Thermocouples  7, 12, and  17  are  approximately  1"  (2.5  cm)  from  adjacent  heaters. 

CROSS  SECTIONAL  VIEW 

0 

INSULATION  LAYER  -  0.75  INCHES  (1.9  CM)  THICK 

_ 1 _ 

INSULATION  LAYER  -  0.75  INCHES  (1.9  CM)  THICK 

_ 2 _ 

INSULATION  LAYER  -  0.75  INCHES  (1 .9  CM)  THICK 

_ 14  13 _ 12 _ 11  10 _ 

_ FLOW-  THROUGH  MODULE  (SEM-E  FORMAT  SIZE) _ 

Picfure  6.  Heater  and  thermocouple  locations. 
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III.  EXPERIMENTAL  PROCEDURES 


A.  DATA  COLLECTION 

Prior  to  conducting  any  experimental  runs,  the  data 
collection  program  was  run  to  determine  the  approximate 
ambient  room  temperature.  In  order  to  prepare  for  data 
collection,  the  refrigerated  bath  was  started  and  the  pump  was 
adjusted  to  a  coolant  flow  level  of  1000  ml/min.  Bath 
temperature  was  adjusted  to  the  ambient  temperature.  The 
power  supply  was  then  energized  and  tuned  to  the  appropriate 
setting  by  running  the  data  collection  program  and  adjusting 
the  voltage  and  current  settings.  Power  settings  are  shown  in 
Table  1.  Because  of  variations  in  resistances  for  each 
heater-precision  resistor  loop,  each  heater  varied  in  its 
power  dissipation.  Thus  the  power  setting  was  based  on  the 
total  power  (Ptotai)  input.  Heater  flux  was  determined  as: 


q 


// 


^total  \ 

_6 _ [ 

1.94 


(3) 


where  1.94  cm^  is  the  effective  heater  area. 

Temperature  data  was  monitored  to  determine  when  steady 
state  had  been  achieved.  Steady  state  was  assumed  when 
thermocouple  readings  on  the  surface  and  above  the  module 
fluctuated  less  than  0 . 15°C  over  a  one  minute  interval .  Once 
at  steady  state,  data  was  recorded  and  saved.  Data  recorded 
by  the  system  included  thermocouple  temperatures  and  voltage 
readings  as  well  as  a  series  of  60  coolant  inlet  and  outlet 
temperature  readings  over  a  period  of  approximately  25 
seconds.  The  flow  rate  was  then  adjusted  to  begin  the  next 
run . 

Once  power  and  flow  rates  were  set,  steady  state  could  be 
reached  in  four  to  thirty  minutes,  depending  upon  the  power 
setting  and  flow  rate.  In  general,  steady  state  was  most 
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rapidly  achieved  at  higher  flow  rates.  Flow  rates  were 
decreased  in  100  ml/min  increments  from  1000  ml/min  to  300 
ml/min.  This  procedure  seemed  to  allow  steady  state  to  be 
reached  more  quickly,  particularly  in  the  insulation  layers. 
Also,  because  the  heaters  were  limited  in  their  ability  to 
withstand  elevated  temperatures  at  high  power  levels,  it  was 
advantageous  to  gradually  raise  surface  temperatures  by 
decreasing  flow,  all  the  while  monitoring  readings  to  avoid 
exceeding  any  design  limitations.  At  flow  rates  lower  than 
300  ml/min  the  flow  was  rather  pulsed,  due  to  the  nature  of 
the  pump,  rather  than  a  steady  flow.  Flow  rates  greater  than 
1000  ml/min  and  less  than  300  ml/min  were  not  used  in  order  to 
avoid  the  extremes  in  flowmeter  settings. 


Table  1.  Power  settings  for  experiments. 

B.  DATA  ANALYSIS 

Data  was  analyzed  as  follows  to  obtain  information  on  the 
heat  dissipation  occurring  at  each  flow  level  and  power 
setting . 
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1.  Calculation  of  Coolant  Thermophysical  Properties 

Thermophysical  properties  were  calculated  by  using  the 
average  of  the  inlet  and  outlet  temperatures  (T(3)  and  T(4) 
respectively)  of  the  coolant  passing  through  the  module, 
determined  as  follows: 


T 


avg 


(  T(4)  +  T(3)  ) 
2 


(4) 


a .  Densi ty  ( g/ml ) 

Manufacturer's  data,  shown  in  the  Table  2,  was  used 
to  determine  the  coolant  density  at  T^^g  by  linear 
interpolation  using  the  equation  shown. 

p  =  0.794  -  (  (8.5E-4)  *  (T^^g  -  20)  )  (5) 


TEMPERATURE 

DENSITY  (g/ml) 

0°C  (32°F) 

0 . 811 

20°C  (68°F) 

0 . 794 

40°C  (104°F) 

0 . 777 

100°C  (212°F) 

0 . 723 

160°C  (320°F) 

0.661 

Table  2.  Brayco  Micronic  889  density  (Castrol,  1994). 


b.  Specific  Heat  Cp  (J/g°C) 

Specific  heat  was  similarly  determined  using  the 
equation  below  and  manufacturer's  data,  shown  in  Table  3, 
through  a  method  of  linear  interpolation. 

Cp  =  (0.52  +  (7.1E-4)  *  (Ta^g  -  10)  )  *  4.184  J/cal  (6) 
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TEMPERATURE 

SPECIFIC  HEAT  (cal/g°C) 

-18°C  (0°F) 

0.49 

10°C  {50°F) 

0.52 

38°C  (100°F) 

0.54 

93°C  (200°F) 

0.58 

Table  3.  Brayco  Micronic  889  specific  heat  (Castrol,  1994) . 


2.  Calculation  of  Heat  Dissipation 

a.  Temperature  Differential 

Sixty  samples  of  inlet  and  outlet  temperatures  (T(3) 
and  T(4) )  taken  at  steady  state  by  the  data  collection  program 
were  used  to  calculate  an  average  delta  T  as  follows: 


60 


AT 


avg 


=  E 


1=1 


Ti{4)  -  Ti(3) 

60 


(7) 


b.  Heat  Dissipated  by  the  Coolant  Qf  (W) 

Q  =  P'^^P^Tavg  (8) 

^  6  0 

where  v  =  flow  rate  in  ml /min 

c.  Heat  Losses  Through  the  Insulation  (W) 

In  order  to  quantify  heat  losses  to  the  atmosphere, 
thermocouple  measurements  in  the  insulation  layers  (see  Figure 
6)  were  used  to  provide  an  estimate  of  conduction  losses  (Qins) 
along  one  heat  flow  path. 
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^  (T(2)  -  T(l)  ) 

JLj 


(9) 


Qir 


where  k 


insulation  conductivity 
0.04  W/  (in°C) , 


A  = 


and  L  = 


surface  area  of  insulation 
5.66  inches  x  6  inches 
33.96  in^,  or  0.0219  , 

insulation  thickness  between  thermocouples 
0.75  inches,  or  0.0191  m. 


The  surface  area  of  the  insulation  was  estimated  as 
the  approximate  area  of  the  recessed  surface  of  the  module 
upon  which  the  heaters  and  thermocouples  were  affixed. 

3.  Calculation  of  Power  into  the  Heaters  P,,  (W) 

Heater  power  was  calculated  by  the  data  collection 
program  using  equations  1  and  2  as  discussed  in  Chapter  II. 
Recorded  values  of  voltage  and  power  were  truncated  to  two 
decimal  places  after  computations  within  the  program  were 
completed. 
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IV.  RESULTS 


A.  RAW  TEMPERATURE  DATA 

Module  surface  and  insulation  layer  temperature  data  at 
the  ten  power  settings  is  shown  in  Tables  4  through  13  (refer 
to  Figure  6  for  thermocouple  location  diagram) .  The  power 
supplied  to  each  heater  and  total  power  into  the  heaters  is 
also  tabulated.  Tables  14  through  23  show  average  inlet  and 
outlet  temperatures,  AT^^g  for  the  coolant,  as  well  as  various 
temperature  differentials  on  the  module  surface.  Graphical 
representations  of  temperature  patterns  on  the  module  surface 
are  shown  in  Figures  7  through  16  for  flow  settings  of  300, 
600,  and  1000  ml /min.  The  average  of  the  inlet  and  outlet 
temperature  differentials  are  plotted  in  Figure  17.  Figure  18 
shows  the  same  data  on  a  log-log  scale. 

B.  HEAT  TRANSFER  DATA 

The  temperature  data  was  input  into  the  equations  in 
Chapter  III  to  produce  the  heat  transfer  calculations  shown  in 
Tables  14  through  23 .  Graphical  representations  of  the  data 
for  Qf  over  the  ten  power  settings  are  shown  in  Figure  19 . 
Independent  laboratory  tests  showed  the  specific  heat  of  the 
coolant  to  be  less  than  that  indicated  by  the  manufacturer  at 
66°C  and  93°C  (based  upon  test  results  received  from  Delsen 
Testing  Laboratories,  Inc.,  on  December  6,  1994).  Figure  20 
shows  a  sampling  of  revised  Qf  calculations  using  an 
extrapolation  of  these  Cp  values.  Chapter  V  discusses  the 
laboratory  information  in  greater  detail,  while  Appendix  B 
includes  a  comparison  of  laboratory  testing  results  and 
manufacturer's  data  for  Cp. 

Figures  21  through  30  illustrate  the  temperature  rise 
experienced  from  "inlet"  to  "exit"  of  each  heater,  (based  on 
the  fluid  flow  path  as  shown  in  Figure  6) ,  for  the  ten  power 
settings  and  a  sampling  of  -flows.  Figures  31  through  40  are 
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graphs  of  the  mid- channel  temperature  drop  experienced  between 
the  "exit"  of  each  heater  and  the  next  thermocouple  placed  on 
the  flow  path  (as  illustrated  in  Figure  6)  .  As  before,  all 
power  levels  are  represented  by  a  sampling  of  flow  settings. 
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HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


RESISTOR  A 


RESISTOR  B 


RESISTOR  C 


RESISTOR  D 


RESISTOR  E 


RESISTOR  F 


300 

400 

8.13 

8.13 

8.12 

8.13 

8.12 

8.12 

8.05 

8.06 

8.12 

8.11 

lEm 

7.98 

0.18 

0.18 

0.18 

0.18 

0.18 

0.17 

0.17 

0.18 

0.18 

0.17 

0.17 

0.17 

0.17 

0.17 

0.18 

0.18 

0.18 

0.17 

0.17 

0.17 

HEATER  A 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

ISSil 

HEATER  B 

2.94 

2.94 

2.94 

2.94 

2.94 

HEATER  C 

2.92 

2.92 

2.92 

2.92 

2.92 

2.92 

2.92 

HEATER  D 

2.86 

2.86 

2.86 

2.86 

2.86 

2.86 

2.86 

HEATER  E 

2.92 

2.91 

2.92 

2.91 

2.91 

2.91 

2.91 

HEATER  F 

2.84 

2.84 

2.85 

2.84 

2.85 

2.85 

2.85 

Total  Power 

17.43 

17.42 

17.44 

17.43 

16.88 

16.92 

17.42 

17.31 

18.04 

17.82 

Thermocouple  temps  (degrees  Celsius 


18.30 


17.73  17.44  17.35  17.33  17.28  17.55 


17.64  17.52  17.48  17.47  1 17.70 


17.37  1 17.06  16.96  16.96  16.94 


18.13  17.78  17.65  17.63  17.56  17.80 


18.27  17.91  17.77  17.73  I  17.65  I  17.90 


17.63 


18.55  1 18.14  1 17.97  1 17.90  17.84  17.99 


17.33  17.28  17.24  17.40 


18.56  18.16 


17.83  17.96 


18.00  I  17.95  I  18.08 


18.59  18.20  18.00 


19.08  18.70  18.48 


18.22  17.83  17.61 


18.27 

18.16 

18.33 

18.22 

Table  4.  Power  and  temperature  data  at  power  setting  1 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


RESISTOR  A 


RESISTOR  B 


RESISTOR  C 


RESISTOR  D 


RESISTOR  E 


RESISTOR  F 


HEATER  A 


HEATER  B 


HEATER  C 


9.37 

9.37 

9.37 

9.37 

9.36 

9.36 

9.37 

9.32 

9.36 

9.35 

9.36 

9.27 

9.27 

9.28 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

iMi 

0.20 

HEATER  E 

3.87 

HEATER  F 

3.84 

3.92 

3.92 

3.92 

3.90 

3.90 

3.90 

3.87 

3.87 

3.88 

3.83 

3.87 

3.87 

3.87 

Total  Power 


18.70 

18.57 

19.24 

19.05 

19.32  I  19.08  18.91 


19.23  I  19.02  I  t9.00  18.89  18.72  18.67 


19.04  18.87  18.83 


18.47  I  18.46  1 18.30 


19.70  I  19.46  I  19.36 


Thermocouple  temps  (degrees  Celsius 


18.13  18.07 


18,64 

18.44 

17.78 

17.69 

19.51  19.31 


19.37  1 19.18  I  18.97 


19.30  19.22 


19.33  19.24 


19.71  1 19.48 


19.76  19.59  1 19.39  I  19.27 


21.02  120.61  120.34  120.12 


19.93 


20.98  20.53  20.22  19.97  19.77  19.65 


21.11  20.62  120.31  120.06  1 19.85  119.73 


Table  5.  Power  and  temperature  data  at  power  setting  2 
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HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


RESISTOR  A 


RESISTOR  B 


RESISTOR  C 


RESISTOR  D 


RESISTOR  E 


RESISTOR  F 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


Total  power 


10.38  10.38 


0.23  0.23 


0.23 


0.23 


0.22 


0.23 


0.23 


Power 


4.89 

4.89 

4.87  : 

4.87 

4.84  1 

4.84 

4.78 

4.78 

4.83 

4.83 

4.80 

4.80 

29.01 

29.01 

.89 

4 

.89 

Thermocouple  temps  (degrees  Celsius) 


20.20  20.23  20.23  20.38 


20.81  120.76  120.77  120.70  120.76 


19.77  I  19.84  19.76  19.58  19.58  1 19.51 


21.72  21.42  21.13  20.96  20.75 


20.57 


20.99  120.88  120.78 


21.07 


21.60  I  21.34  121.22 


20.96 


21.08 


11 

23.24 

21.62 

21.44  21.29 

21.29 

12 

22.09 

IBfEIBBII 

20.87 

20.61 

20.49  20.34 

20.35 

22.68  22.24  21.94  21.65 


22.71  22.24 


23.49  22.99  22.64  22.36  22.20 


21.32  21.01  20.86 


22.21  22.01 


123.60  123.01  122.64  122.31  122.13 


21.47 


22.00 


21.44 


21.31 


21.93 


Table  6.  Power  and  tempe-rature  data  at  power  setting  3 
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Flow  settina  (ml/min 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


RESISTOR  A 


RESISTOR  B 


RESISTOR  C 


RESISTOR  D 


RESISTOR  E 


RESISTOR  F 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


300 

400 

500 

600 

700 

800 

|g!a.lKM.M 

Voltages  (V) 

11.51 

11.51 

11.51 

11.51 

11.51 

11.51 

11.51 

11.51 

UiUIBI.-CIRIifilKBI.-fil 

11.50 

11.51 

11.50 

11.50 

11.40 

11.40 

11.49 

11.49 

11.49 

11.49 

11.49 

UnBslIEBUlICgtgl 

11.30 

11.30 

11.28 

11.28 

11.29 

11.28 

11.26 

11.27 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

E^l 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.24 

0.24 

0.24 

0.24 

0.24 

0.24 

EEBl 

0.24 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

itISI 

0.25 

0.25 

0.25 

0.24 

0.25 

0.25 

0.24 

0.24 

0.24 

Power  (W) 

5.90 

5.90 

5.90 

5.90 

5.90 

5.90 

!  5.90 

5.90 

5.87 

5.87 

5.87 

5.88 

5.88 

5.88 

EKSi 

5.87 

5.83 

ESS 

5.83 

5.84 

ESI 

5.84 

5.84 

5.83 

5.71 

5.72 

5.72 

5.72 

5.72 

5.72 

5.72 

5.72 

5.83 

5.83 

5.83 

5.83 

5.83 

5.84 

5.83 

5.83 

5.68 

5.68 

WttM 

5.67 

5.68 

5.67 

5.65 

5.65 

34.82 

34.84 

34.82 

ThermocouDle  temos  (deorees  Celsius 


19.54  19.53 


20.60  20.48  20.44  20.37  20.32 


18.52  I  18.67  18.69  18.82  18.92 


21.56  20.98  20.68  20.66  20.53  20.32 


20.38  20.30  20.24  20.30  20.32  20.15 


20.58 

20.45 

20.74 

20.61 

20.52 

20.36 

Table  7.  Power  and  temperature  data  at  power  setting  4 
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Flow  settinq  (ml/min 


Heater  A 


Heater  B 


Heater  C 


Heater  D 


Heater  E 


Heater  F 


RESISTOR  A 


RESISTOR  B 


300 

400 

500  600 

Voltages  (V) 

^^1 

12.41 

12.41  12.41 

12.41 

12.41 

12.41  12.40 

12.40 

12.40 

12.40 

12.40  12.40 

12.40 

12.40  12.40 


12.34  I  12.34 


12.39  12.39  12.39  12.39  12.39  12.39 


12.30  12.30  12.30  12.30  12.30 


RESISTOR  C 

0.27 

0.27 

0.27 

0.27 

0.27 

0.27 

0.27 

0.27 

RESISTOR  D 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 

RESISTOR  E 

0.27 

0.27 

0.27 

ItBM 

0.27 

0.27 

0.27 

0.27 

RESISTOR  F 

0.27 

0.27 

0.27 

0.27 

WtSM 

0.27 

0.27 

0.27 

Power  (W) 


Heater  A 

6.85 

6.85 

6,85 

6.85 

6.85 

6.85 

Heater  B 

6.82 

6.82 

6.82 

6.82 

6.82 

6.82 

6.82 

6.82 

Heater  C 

6.78 

6.78 

6.78 

6.78 

6.78 

6.78 

6.78 

6.78 

Heater  D 

6.70 

6.70 

6.70 

6.70 

6.70 

6.70 

6.70 

6.70 

Heater  E 

6.77 

6.77 

6.77 

6.77 

6.77 

6.77 

6.77 

6.77 

Heater  F 

6.73 

6.73 

6.73 

6.73 

6.73 

6.73 

Total  Power 

1  40.65 

40.65 

40.65 

40.65 

40.65 

40.65 

21.35  21.19 

21.08 

22.99  22.76 

22.63 

25.08  24.34  23.89  23.50 


Table  8.  Power  and  temperature  data  at  power  setting  5 
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HEATER  A 

13.29 

HEATER  B 

13.28 

HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


13.29 

13.29 

13.29 

13.28 

13.28 

13.28 

13.28 

13.28 

13.27 

13.28 

13.26  13.26 


13.07  13.07 


13.26  13.26 


13.07  13.09 


RESISTOR  A 

0.29 

0.29 

0.29 

0.29 

RESISTOR  B 

0.29 

0.29 

0.29 

0.29 

0.29 

0.29 

0.29 

0.29 

RESISTOR  C 

0.29 

0.29 

0.29 

0.29 

0.29 

0.29 

0.29 

Itl^ 

RESISTOR  D 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

RESISTOR  E 

0.29 

0.29 

0.29 

0.29 

0.29 

l»K3i 

0.29 

0.29 

RESISTOR  F 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

\mm^ 

7.82 

7.82 

7.76 

7.77 

7.62 

7.62 

Power 


7.86  I  7.86 


ThermocouDle  temps  {deqrees  Celsius 


18.69 


19.46  19.47  19.50  19.38 


21.04  20.79  20.65  20.43  20.16  19.95 


18.30  18.34  18.22  18.41  18.42  I  18.44 


20.39  20.21 


20.25  20.15  20.15 


21.19 


19.63  19.34  19.11 


21.66  21.24  20.92 


22.04  21.58  21.24  21.16  21.01  20.91 


22.05  21.45 


20.56 


22.81  22.18  121.74  121.52  121.36  121.19 


21.03  20.45 


22.96  22.27  21.79  21.57  121.37  121.21 


22.61 


18.81 


20.41 


20.64 


20.31 


20.94 


23.69 

21.34 

23.35 

26.29 

24.47  1 

23.60 

Table  9.  Power  and  tempexature  data  at  power  setting  6 
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400 

500  600 

700  i 

800  i 

900 

1000 

HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


Total  Power 


Voltages  (V) 


14.08  14.08  14.08  14.08  14.08  14.08 


14.08  14.08  14.07  14.07  14.07  14.08 


14.07  14.07  14.07  14.07  14.07  14.07 


14.00  1 13.99  14.00  1 14.00  14.00 


14.06  14.06  14.06  14.05  14.05  1 14.06 


13.93  13.94  13.93  13.94  13.94  13.94 


RESISTOR  A 

0.31 

0,31 

0.31 

0.31 

EEOIEEfli 

0.31 

0.31 

RESISTOR  B 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

RESISTOR  C 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

RESISTOR  D 

0.30 

0.30 

0.30 

0.30 

EISI 

0.30 

RESISTOR  E 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

RESISTOR  F 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

Power  (W) 


8.81 


8.77 


8.71 


8.61 


8.71 


8.63  8.64 


52.25  I  52.25 


Thermocouple  temps  (degrees  Celsius 


20.33  119.20 


21.12 


21.94  21.68  21.32 


18.92  18.98  18.95  19.02  19.01 


22.46  21.91  21.62  121.31  121.02 


21.82  21.59  21.54  21.39  121.13 


22.13  21.94  21.69 


21.82  21.62  21.34 


22.58  22.35  22.08 


20.77  20.57  120.39 


22.65  22.42 


23.00  1  22.72 


22.72  22.46 


25.63  1 

24.80 

24.28 

23.93  i 

23.65  23.42  23.01 

24.09  23.63  23.35  23.15  22.69 


25.88  124.90  124.31  23.85  23.55  123.34  122.86 


Table  10.  Power  and  temperature  data  at  power  setting  7, 
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Flow  settinqs  (ml/min 


HEATER  A 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

HEATER  B 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

14.84 

HEATER  C 

14.84 

14.84 

14.84 

T4.83 

14.83 

14.83 

14.83 

14.84 

HEATER  E 


HEATER  F 


RESISTOR  A 


RESISTOR  B 


RESISTOR  C 


RESISTOR  D 


14.72  14.73 


0.32  0.32 


14.71 

14.70 

0.32 

0.32 

0.32 

0.32 

0.32 

0.32 

0.32 

0.32 

mmi 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


Total  Power 


Power  (W 


9.79 

9.79 

9.79 

9.79 

9.75 

9.75 

9.75 

9.75 

9.68 

9.68 

9.68 

9.68 

9.68 

9.68 

9.54 

9.54 

9.54 

9.54 

9.54 

9.54 

9.68 

9.68 

9.68 

9.68 

9.68 

9.68 

9.62 

9.62 

9.62 

9.62 

58.06 

58.07 

58.06 

58.06 

19.01 


20.11 


21.34 


19.23 


21.38 


Thermocouple  temps  (degrees  Celsius) 


18.79  18.69  18.68  18.71 


19.96 


21.41  121.36 


19.17  19.26  i  19.33 


22.11 


21.36  21.34  21.31 


24.16  23.56  23.17 


27.07  25.43  24.57  23.95  23.56 


26.65  25.05  24.20  23.58  123.22 


28.29  26.57  25.64  24.99  24.56 


22.97  22.34  21.94 


24.36 


28.92  26.91  25.82  25.07  24.56 


24.08 


21.49 


23.81 


23.99 


22.85 


22.58 


23.81 


21.27 


23.50 


23.69 


Table  11.  Power  and  temperature  data  at  power  setting 
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HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


RESISTOR  A 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

RESISTOR  B 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

RESISTOR  C 

0.33 

0.33 

0.33 

0.33 

0.33 

RESISTOR  D 

0.33 

ItBcB 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

RESISTOR  E 

0,33 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

RESISTOR  F 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

HEATER  D 


HEATER  E 


HEATER  F 


Total  Power 


Power  (W) 


HEATER  A 

KBKHlIBICTiHCT 

10.77 

10.77 

10.77 

10.77 

10.77 

HEATER  B 

10.72 

10.72 

10.72 

10.72 

10.72 

HEATER  C 

10.65 

10.65 

10.65 

17.74 

18.11 

18.32  1 

WSIBm 

18.28 

18.26 

18.33 

25.87 

22.27 

21.54 

21.21 

23.06 

22.22 

21.88 

20.94 

20.10 

19.73 

24.02 

22.90 

22.39 

24.68 

23.44 

22.83 

25.00 

23.40 

22.65 

26.13 

24.45 

23.65 

26.59 

24.76 

23.84 

28.38 

26.39 

25.41 

25.35 

29.16 

26.77 

25.62 

Thermocouple  temps  (decrees  C 


18.23  18.17 


19.60  I  19.41  19.32 


20.89 


18.39  18.23  18.34 


21.98  21.53  21.17 


21.02  20.69  I  20.59 


21.59 


19.54  I  19.27 


22.05  21.72  21.59 


21.93 


22.40 


20.85  120.47  20.27 


22.73  22.46 


23.65  123.16  22.86 


21.87  21.40 


19.23 

19.18 

20.81 

20.84 

18.11 

18.25 

18.93 

18,99 

21.19 

21.18 

21.50 

21.50 

21.07 

20.98 

21.95 

21.88 

19.84 

19.78 

21.97 

21.85 

22.37 

22.25 

22.04 

21.89 

23.41 

23.24 

23.13 

22.94 

23.35 

23.14 

Table  12 .  Power  and  temperature  data  at  power  setting  9 
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26.05  25.35  24.50  24.29  23.93  23.75 


26.29  24.07  22.91  22.23  21.39  21.23 


29.80  27.28  25.99  25.17  24.29  24.00  123.62  123.39 


30.19  127.62  126.26  125.45  124.51  124.27  123.85  23.64 


Table  13.  Power  and  temperature  data  at  power  setting  10 
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DEGREES  CELSIUS 


Figure  7.  Module  surface  'temperatures  at  power  setting  1. 


Figure  8.  Module  surface .temperatures  at  power  setting  2. 
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DEGREES  CELSIUS 


Figure  9.  Module  surface  temperatures  at  power  setting  3. 
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DEGREES  CELSIUS 


Figure  10.  Module  surface  temperatures  at  power  setting  4. 
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DEGREES  CELSIUS 


Figure  11.  Module  surface- temperatures  at  power  setting  5. 
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Figure  12.  Module  surface 


mperatures  at  power  setting  6 . 


2 


DEGREES  CELSIUS 


DEGREE 


THERMOCOUPLE  NUMBER 


Figure  14.  Module  surface  temperatures  at  power  setting  8. 
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DEGREES  CESLIUS 


Figure  15.  Module  surface- temperatures  at  power  setting  9. 
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DEGREES  CELSIUS 


Figure  17.  Coolant  temperature  increase  through  FTM. 
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10 


■  POWER  SETTING  1 

❖  POWER  SETTING  2 
£  POWER  SETTING  3 
□  POWER  SETTING  4 
e  POWER  SETTINGS 
A  POWER  SETTINGS 

•  POWER  SETTING  7 

★  POWER  SETTING  8 
▼  POWER  SETTING  9 
©POWER  SETTING  10 


1000 

ml/min 


temperature  differentials. 


T(3)avg  (inlet) 


T(4)avg  (outlet) 


delta  Tfluid  (avg) 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


16.81 

16.81 

16.77 

16.75 

16.68 

16.61 

16.98 

18.89 

18.34 

17.95 

17.73 

17.59 

mggimrViii 

2.08 

1.53 

0.90 

0.83 

0.72 

Delta  T  Across  Heaters  (degrees  C) 

0.26 

0.23 

0.20 

0.17 

0.15 

0.19 

0.14 

0.17 

0.14 

0.13 

0.12 

0.10 

0.09 

0.11 

0.30 

0.29 

0.26 

0.26 

0.27 

0.26 

0.26 

0.16 

0.16 

0.11 

0.12 

0.10 

0.13 

0.11 

0.58 

0.49 

0.50 

0.48 

0.45 

0.44 

0.44 

0.08 

0.05 

0.05 

0.06 

0.06 

Mid-channel  Temperature  Drops  (degrees  C) 


After  heater  A 

0.60 

0.59 

0.58 

0.56 

0.52 

0.54 

0.53 

lyyw 

After  heater  C 

0.69 

0.68 

0.64 

0.64 

0.62 

0.60 

0.58 

0.56 

After  heater  E 

0.84 

0.86 

0.86 

0.87 

0.86  1 

0.83 

0.83 

0.85 

Coolant  Thermophysical  Properties 


density  (g/ml) 


Cp  (J/(gC)) 


mm 

0.7964 

0.7965 

2.197 

2.197 

2.197 

Qf  (W) 

18.20 

17.85 

17.21 

Qins  (W) 

0.03 

0.02 

0.02 

q"ins  (W/m) 

1.31 

1.06 

0.88 

Heat  Dissipation 


Table  14.  Calculations  for  power  setting  1. 


49 


0.15 

0.15 

0.14 

0.12 

0.10 

0.11 

0.33 

0.33 

0.31 

0.15 

0.14 

0.17 

0.54 

0.54 

0.53 

0.08 

0.08 

0.10 

Mid-channel  Temperature  Drops  (degrees  C) 


After  heater  A 

0.77 

0.72 

0.75 

0.73 

0.71 

0.67 

After  heater  C 

0.93 

0.87 

0.85 

0.84 

0.80 

0.77 

0.78 

0.74 

After  heater  E 

1.11 

1.10 

1.12 

1.09 

1.07 

1.07 

1.06 

1.05 

Coolant  Thermophysical  properties 


density(g/mi) 

0.7954 

0.7954 

0.7955 

Cp  (J/(gC)) 

2.202 

2.201 

2.201 

Heat  Dissipation 


Qf{W) 


Qins  (W) 


q”ins  (W/m) 


23.81 

23.69 

24.04 

24.16 

25.38 

0.02 

0.02 

0.02 

0.02 

0.02 

1.00 

0.96 

0.96 

0.88 

0.71 

Table  15.  Calculations  for  power  setting  2. 
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T{3)avg  (inlet) 


T(4)  avg  (outlet) 


deltaTfluid  (avg) 


HEATER  A 

0.42 

0.36 

HEATER  B 

0.27 

0.23 

21.43  21.08  20.89  20.76 


1.74  1.47  1.35  1.18 


2.51  2.01 


Delta  T  Across  Heaters  (degrees  C) 


19.69 


20.77 


HEATER  D 


HEATER  E 


HEATER  F 


0.26 

0.20 

0.16 

0.16 

0.44 

0.39 

0.18 

0.19 

0.66 

0.68 

0.10 

0.12 

Mid-channel  Temperature  Drops  (degrees  C) 


After  heater  A 

1.00 

1.02 

0.98 

0.94 

0.92 

0.85 

After  heater  C 

1.15 

1.09 

1.06 

1.02 

1.01 

0.95 

0.95 

0.94 

After  heater  E 

1.38 

1.38 

1.36 

1.32 

1.35 

1.34 

1.33 

1.30 

density  (g/ml) 


Cp  (J/(gC)) 


Qf(W) 


Qins  (W) 


q"ins  (W/m) 


Coolant  Theimophysical  Properties 


2210  1  2.208  12.208  1  2.207  1 2207  2.206  2.206  2.206 


28.91  29.31 


0.05  0.03 


30.47 

30.04 

31.52 

0.03 

0.02 

0.02 

1.16 

0.97 

0.93 

Table  16.  Calculations  for  power  setting  3. 
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2.96  J  2.33  2.06  1.71  1.61  |  1.41 


Delta  T  Across  Heaters  (degrees  C) 


HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


0.29 

023 

0.59 

0.57 

0.27 

0.25 

1.02 

0.96 

0.20 

0.17 

Mid-channel  Temperature  Drops  (degrees  C) 


1.19 

1.18 

1.15 

1.13 

1.10 

1.11 

1.03 

1.38 

1.32 

1.29 

1.24 

1.18 

1.16 

1.14 

.75 

1.73 

1.72 

1.70 

1.67 

1.67 

1.65 

density  (g/ml) 


Cp  (J/(gC)) 


Qf(W) 


Qins  (W) 


q"lns  (W/m) 


0.7935 


2.207 


_ Coolant  Thermophysical  Properties 


0.7940  0.7941  0.7943  0.7942  0.7943  0.7944 


12.206  2.205  2.205  2.205 


_  Heat  Dissipation 


34.00  36.08  34.94  37.60  37.03 


0.05  0.04  0.04 


1.91 


36.77 


0.03 


1.56 


Table  17.  Calculations  for  power  setting  4 
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Table  18.  Calculations  for  power  setting  5. 
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Table  19.  Calculations  for  power  setting  6. 
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Flow  setting  (ml/min) 


300  400  500  600  700  800  900  1000 


Total  Power  Input  (W) 


52.26  52.26  52.27  52.24  52.25  52.25  52.28  52.26 


Degrees  C 


T(3)avg  (inlet) 

18.54 

18.85 

18.93 

18.94 

19.04 

18.97 

18.77 

18.66 

T(4)avg  (outlet) 

delta  Tfluid  (avg) 

6.07 

4.45 

3.54 

3.02 

2.64 

2.37 

2.16 

1.97 

Delta  T  Across  Heaters  (degrees 

C) 

HEATER  A 

0.55 

0.47 

0.38 

0.36 

eesesi 

0.34 

0.28 

HEATER  B 

0.51 

0.42 

0.36 

0.33 

0.30 

0.28 

0.27 

0.26 

HEATER  C 

0.92 

0.78 

0.76 

0.73 

0.74 

0.71 

HEATER  D 

0.39 

0.34 

0.35 

l»WmcgH{|cf| 

HEATER  E 

1.53 

mMwm 

1.31 

1.20 

1.19 

1.22 

1.12 

HEATER  F 

0.31 

0.26 

0.22 

0.22 

0.22 

0.20 

0.19 

0.17 

density  (g/ml) 


Cp  (J/(gC)) 


Of  (W) 


Qins  (W) 


q"ins  (W/m) 


Mid-channel  Temperature  Drops  (degrees  C) 


After  heater  A 

1.71 

.1.69 

1.68 

1.61 

After  heater  C 

2.05 

1.99 

1.91 

1.86 

1.81 

After  heater  E 

2.44 

2.45 

2.46 

2.43 

2.39 

2.35 

Coolant  Thermophysical  Properties 


0.7927  0.7931  0.7934  0.7936  0.7937  |0.7939  |0.7941  0.7943 


2.208  2.207  2.207  2.206  2.205  2.205 


Heat  Dissipation 


53.19 

51.97 

51.68 

0.10 

0.08 

0.08 

4.72 

3.80 

3.44 

Table  20.  Calculations  for  power  setting  7. 
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HEATER  D 


HEATER  E 


HEATER  F 


Mid-channel  Temperature  Drops  (degrees  C) 


After  heater  A 

1.93 

1.95 

1.93 

1.87 

1.84 

1.81 

1.76 

1.71 

After  heater  C 

2.27 

2.21 

2.13 

2.06 

1.99 

1.96 

1.92 

1.90 

After  heater  E 

2.64 

2.67 

2.66 

2.66 

2.61 

2.59 

2.60 

2.54 

Qf{W) 


Qins  (W) 


q"ins  (W/m) 


0.7918 


2.213 


57.87 


0.11 


4.89 


Coolant  Thermophysicai  Properties 


0.7925  0.7929  0.7932  0.7934  0.7935  0.7937 


2.209  2.208  1 2.208  2.207  2.207 


Heat  Dissipation 


56.61  57.16  58.30  58.72  59.99 


0.7938 


2.206 


Table  21.  Calculations  for  power  setting  8. 


56 


T(3)  avg  (inlet) 

18.31 

18.33 

T(4)avg  (outlet) 

25.84 

23.75 

deltaTfluid(avg) 

7.53 

5.43 

HEATER  A 


HEATER  B 


HEATER  C 


HEATER  D 


HEATER  E 


HEATER  F 


0.68 

0.67 

0.57 

0.60 

0.63 

0.54 

0.44 

0.43 

0.38 

0.34 

1.05 

1.00 

■  0.96 

0.94 

0.92 

0.49 

0.46 

0.45 

0.42 

0.40 

1.63 

1.57 

1.51 

1.45 

1.44 

0.32 

0.27 

0.27 

0.20 

0.23 

Mid-Channel  Temperature  Drops  (degrees  C) 


After  heater  A 

2.13 

2.13 

2.14 

2.06 

2.02 

1.99 

1.92 

1.88 

After  heater  C 

2.51 

2.42 

2.33 

2.25 

2.18 

2.13 

2.11 

2.10 

After  heater  E 

2.83 

2.89 

2.87 

2.87 

2.83 

2.80 

2.77 

2.72 

density  (g/ml) 


Cp  (J/(gC)) 


Qf  (W) 


Qins  (W) 


q"ins  (W/m) 


0.7922 


2.212 


65.99 


Coolant  Thermophysical  Properties 


0.7931  0.7935  0.7938  0.7941  0.7943  0.7945  0.7945 


12.209  2.207  2.206  1 2.205  2.205 


Heat  Dissipation  . 


66.86 

67.72 

67.43 

68.90 

0.08 

0.07 

0.07 

0.08 

3.41 

3.30 

3.31 

3.48 

Table  22 


Calculations  for  power  setting  9 . 
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After  heater  A 

2.33 

2.33 

2.30 

2.24 

2.18 

2.19 

2.10 

waam 

After  heater  C 

2.72 

2.64 

2.57 

2.49 

2.39 

2.33 

2.30 

After  heater  E 

3.12 

3.14 

3.15 

3.12 

3.12 

3.06 

3.01 

IBS 

Coolant  Thermophysical  Properties 


density  (g/ml) 

0.7918 

0.7928 

0.7934 

0.7937 

0.7942 

0.7942 

0.7944 

0.7945 

Cp  (J/(gC)) 

2.208 

2.207 

2.205 

2.205 

2.204 

2.204 

Heat  Dissipation 


Qf(W) 

71.14 

69.15 

69.64 

71.12 

71.71 

72.85 

Qins  (W) 

0.12 

0.10 

0.09 

0.08 

0.06 

0.08 

0.07 

0.07 

q"ins  (W/m) 

5.39 

4.61 

4.24 

3.82 

2.50 

3.49 

3.39 

2.97 

Table  23.  Calculations  for  power  setting  10. 


Figure  19 .  Heat  absorbed  by  fluid  (Qj )  . 


Qf(W) 


Figure  20.  Revised  Qf  calculations  based  on  lab  data, 
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DEGREES  CELSIUS 


Figure  21.  Temp  rise  across  heaters  at  power  setting  1. 
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DEGREES  CELSIUS 


Figure  23.  Temp  rise  across  heaters  at  power  setting  3. 


63 


CELSIUS 


Figure  24.  Temp  x -  ^e  across  heaters  at  power  setting  4. 
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DEGREES  CELSIUS 


POWER  SETTING  5 


Figure  25.  Temp  rise  across  heaters  at  power  setting  5. 
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Figure  26.  Temp  rise  across  heaters  at  power  setting  6 
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DEGREES  CELSIUS 


POWER  SETTING  7 


Figure  27.  Temp  rise  across  heaters  at  power  setting  7. 
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DEGREES  CELSIUS 


Figure  28.  Temp  r.,ae  across  heaters  at  power  setting  8. 
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DEGREES  CELSIUS 


POWER  SETTING  9|| 
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Figure  29.  Temp  rise  across  heaters  at  power  setting  9. 
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DEGREES  CELSIUS 


POWER  SETTING  10 


Figure  30.  Temp  rise  across  heaters  at  power  setting  10. 
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DEGREES  CELSIUS 


Figure  31.  Mid- channel -temp  drops  at  power  setting  1 
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DEGREES  CELSIUS 


POWER  SETTING  2 


iSOOML/MIN 
D  600  ML/MIN 
1000  MUMIN 


After  heat 


After  heater  > 


After  heater  E 


Figure  32. 


M..  ..-channel  -temp  drops  at  power  setting  2. 
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DEGREES  CELSIUS 


Figure  33.  Mid-channel  -temp  drops  at  power  setting  3 


DEGREES  CELSIUS 


Figure  34.  Mid-channel  -temp  drops  at  power  setting  4. 
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Figure  35.  Mid-channel  -temp  drops  at  power  setting  5 
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DEGREES  CELSIUS 


DEGREES  CELSIUS 


3 


2.5 


1.5 


1 


S  300  MUMIN 
n  600  ML/MIN 
■  1000  ML/MIN 


Figure  37.  Mid-channel -temp  drops  at  power  setting  7. 


DEGREES  CELSIUS 


Figure  38.  Mid-channel  -temp  drops  at  power  setting  8. 
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DEGREES  CELSIUS 


Figure  39.  Mid-channel  -temp  drops  at  power  setting  9 
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DEGREES  CELSIUS 


Figure  40.  Mid-channel  temp  drops  at  power  setting  10. 
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V.  DISCUSSION  AND  CONCLUSIONS 


A.  TEMPERATURE  PATTERNS 

The  data  presented  in  Tables  4  through  13  as  well  as  the 
graph  in  Figure  17  most  clearly  portray  the  cooling 
capabilities  of  the  module.  As  might  be  expected  after 
viewing  Figure  3  (illustrating  the  coolant  flow  path)  and 
Figure  6  (showing  thermocouple  locations) ,  temperatures  tend 
to  rise  across  each  heater  and  drop  in  the  flow  channel 
between  heaters.  Figures  21  through  40  show  these  tendencies. 
The  temperatures  on  the  FTM  surface  basically  exhibit  an 
increasing  nature  along  the  flow  path,  with  the  exception  of 
data  from  the  three  thermocouples  located  between  heaters 
(T(7) ,  T(12) ,  and  T(17) ) . 

It  also  appears  that  as  the  coolant  enters  the  module,  it 
has  a  cooling  effect  on  the  flow  immediately  to  its  right  as 
viewed  in  Figure  6.  The  exact  nature  of  the  passage  as  the 
fluid  turns  from  one  channel  to  the  next  is  unclear,  but  one 
would  suspect  behavior  in  this  area  to  be  similar  to  that  of 
a  counterflow  heat  exchanger  due  to  the  opposing  directions  of 
the  flow  inside  the  module  and  in  the  entrance  channel.  The 
results  show  a  pattern  of  relatively  smaller  temperature 
increases  across  heater  D  as  measured  by  thermocouples  13  and 
14,  as  well  as  a  slight  decrease  at  thermocouple  15,  where  the 
flow  makes  the  turn  into  the  final  flow  path  through  the  FTM. 
Although  the  flow  makes  a  similar  turn  in  proceeding  from 
thermocouple  9  to  thermocouple  10,  the  temperature  continues 
to  increase  in  that  area  at  the  lower  flow  rates. 

The  lowest  temperature  measured  on  the  module  surface 
invariably  occurred  at  thermocouple  7 ,  located  between  heaters 
in  the  first  flow  pass.  The  highest  surface  temperature 
varied,  but  was  consistently  at  either  thermocouple  16  or  19, 
located  above  the  final  pass  after  heaters  E  and  F 
respectively.  The  temperature  midway  in  this  flow  pass,  as 


81 


measured  by  thermocouple  17,  indicated  a  sharp  decrease  and 
thus  provides  a  strong  testimony  to  the  effectiveness  of  the 
internally  occurring  convection,  at  least  in  this  particular 
area. 

While  efforts  were  made  to  ensure  consistency  and 
accuracy,  the  narrow  range  of  temperatures  being  measured 
created  difficulties.  Slight  flutter  in  the  flow  rate 
combined  with  temperature  variations  exiting  the  reservoir  in 
the  refrigerated  bath  caused  noticeable  fluctuations  in  the 
coolant  inlet  temperature  and  thus  the  surface  temperature 
data.  Such  fluctuations,  which  approached  0.2°C  at  the  higher 
flow  and  the  higher  power  settings,  indicated  that  the  bath 
was  experiencing  difficulty  in  uniformly  cooling  the  flow 
prior  to  its  return  to  the  module.  This  provided  the 
motivation  for  averaging  sixty  of  the  inlet  and  outlet 
temperatures  over  approximately  25  seconds.  Still,  it  must  be 
recognized  that  these  two  temperatures  would  have  to  be 
precisely  spaced  in  order  to  ensure  that  the  measurements 
reflected  both  the  inlet  and  outlet  temperatures  of  the  same 
"slug"  of  fluid.  Because  of  their  relatively  simultaneous 
measurement,  there  is  an  inherent  inaccuracy  in  the  resulting 
temperature  differential. 

It  was  noted  during  the  experiments  that  a  slight 
increase  (or  decrease)  in  inlet  temperature  for  runs  at  the 
same  power  level  and  flow  rate  would  result  in  relatively  the 
same  magnitude  increase  (or  decrease)  in  module  surface 
temperatures,  although  the  temperature  differential 
experienced  between  the  inlet  and  outlet  was  relatively 
stable.  For  example,  decreasing  inlet  temperature  0.5°C  would 
result  in  a  similar  drop  in  surface  and  exit  temperatures,  but 
the  temperature  rise  from  inlet  to  outlet  would  not 
significantly  change.  With  this  knowledge,  efforts  were  made 
to  be  relatively  consistent  regarding  inlet  temperatures.  It 
was  observed  that  if  the  flow  were  entering  at  a  temperature 
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higher  than  room  temperature,  insulation  layer  temperatures 
might  not  steadily  decrease  from  the  inner  layer  (thermocouple 
2)  to  the  outer  layer  (thermocouple  0)  ,  indicating  the 
likelihood  of  heating  from  the  ambient  air.  Such  heating 
could  only  be  assumed  to  skew  all  data.  In  general,  inlet 
temperatures  were  sought  in  the  vicinity  of  room  temperature 
as  determined  by  readings  prior  to  energizing  the  refrigerated 
bath  or  power  supply. 

In  some  instances,  the  data  may  show  surface  temperatures 
to  have  actually  decreased  when  in  fact  power  was  held 
constant  and  flow  was  decreased,  creating  a  situation  where 
the  opposite  would  be  expected.  The  absence  of  a  temperature 
rise  on  the  module  surface  may  be  best  explained  by  noting  the 
inlet  temperature  (thermocouple  3)  and  the  temperature  of 
thermocouple  0,  which  was  near  ambient.  If  an  odd  decrease  in 
surface  temperatures  occurs  in  the  data,  these  two 
thermocouples  tend  to  indicate  that  either  entrance 
temperature  or  room  temperature  had  increased  (or  both) . 

During  the  course  of  experimental  runs,  the  room 
temperature  undoubtedly  changed.  Bath  temperature  was 
frequently  adjusted  as  well  since  at  higher  flow  rates  the 
ability  of  the  bath  to  cool  fluid  prior  to  its  re-entering  the 
flow  was  challenged. 

B.  HEAT  TRANSFER  CALCULATIONS 

Though  the  temperature  differential  between  fluid  inlet 
and  outlet  steadily  increased  with  decreasing  flow  rate  as 
expected  (Figure  17) ,  the  energy  absorbed  by  the  coolant  (Q^) 
did  not  always  decrease  with  decreasing  flow  (Figure  19)  . 
Even  more  unusual  was  the  tendency  to  have  the  calculated 
exceed  the  total  power  to  the  heaters.  This  occurred  in  60  of 
the  80  runs  (75%),  with  calculations  yielding  up  to  a  13.8% 
increase.  Similar  baffling  results  were  produced  in 
experiments  conducted  by  the  Crane  Division  of  the  Naval 
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Surface  Warfare  Center,  where  16  of  23  (69.6%)  flow  and  power 
settings  resulted  in  thermal  efficiencies  in  excess  of  100%, 
with  values  as  high  as  110.9%  calculated  (Buechler  and  Brough, 
1993),  though  no  conclusion  was  reached  as  to  why  this  was 
occurring. 

Several  possible  avenues  for  resolving  these  results  were 
explored.  Because  each  calculation  of  heat  absorbed  by  the 
fluid  depended  upon  only  three  experimentally  measured 
quantities  (flow  and  inlet  and  outlet  temperature) , 
calibration  of  the  flowmeter  and  thermocouples  3  and  4  was 
repeated  to  verify  previous  results.  While  the  thermocouple 
recalibration  using  the  Rosemount  92 OA  commutating  bridge  and 
a  constant  temperature  ethylene  glycol  bath  narrowed  the 
uncertainty,  the  flowmeter  calibration  fell  within  the 
original  uncertainty.  Uncertainty  calculations  (see  Appendix 
B)  took  into  account  flowmeter  calibration  and  temperature 
uncertainty,  however  the  values  obtained  were  not  sufficient 
to  account  for  the  discrepancy  in  all  cases. 

Data  was  collected  over  a  two  hour  period  beginning  when 
steady  state  appeared  to  have  been  reached  for  the  same  power 
and  flow  settings  to  determine  if  inadequate  time  was  being 
allowed  between  runs.  No  noticeable  changes  in  temperature 
data  were  encountered. 

Thermocouple  placement  in  the  inlet  and  exit  flow  was 
also  adjusted  with  no  change  in  results.  The  thermocouples 
were,  however,  approximately  11.5  cm  from  the  actual  inlet  and 
exit  of  the  module  due  to  limitations  in  the  assembly  of 
tubing  and  connections.  This  could  very  possibly  account  for 
a  small  temperature  difference  that  could  bring  the  data  into 
a  feasible  range.  Because  of  the  relatively  small  power 
levels,  slight  temperature  differences  had  significant  impacts 
on  results. 

Coolant  properties  (density  and  specific  heat)  used  in 
calculations  were  based  on  linear  interpolations  between 
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values  supplied  by  the  manufacturer.  A  hydrometer  was  used  to 
measure  coolant  density  at  20°C  as  0.7975  g/ml,  while  the 
manufacturer '^s  data  specified  a  density  of  0.794  g/ml.  While 
possibly  a  very  significant  discovery,  this  would  only  tend  to 
drive  the  measured  heat  dissipation  levels  further  beyond  100% 
thermal  efficiency. 

Having  pursued  several  possible  solutions  to  no  avail, 
the  specific  heat  of  the  fluid  was  investigated  to  determine 
whether  it  matched  manufacturer's  data.  Testing  performed  by 
an  independent  laboratory  using  a  differential  scanning 
calorimeter  gave  a  possible  indication  that  specific  heat  data 
used  for  the  experiments  was  erroneous.  Two  tests  at  93°C 
resulted  in  Cp  values  of  0.57  cal/ (g°C)  ,  slightly  less  than 
manufacturer's  data  shown  in  Table  3.  Two  tests  conducted  at 
66°C  (150°F)  indicated  C  values  of  0.54  cal/ (g°C)  ,  while  a 
linear  interpolation  using  Table  3  produces  0.56  cal/(g°C). 
Assuming  that  the  specific  heat  is  indeed  linear  in  this 
region,  as  Table  3  seems  to  indicate,  an  extrapolation  of  the 
test  data  and  utilization  of  the  obtained  figures  in 
determining  specific  heat  would  lower  calculated  specific  heat 
values  by  approximately  seven  percent  in  the  range  of  testing. 
As  a  result,  energy  balance  calculations  for  the  coolant  would 
also  be  reduced  by  the  same  factor.  Thus,  the  76.17  W  value 
obtained  for  power  setting  10  (69.68  W)  at  1000  ml/min  would 
become  70.88  W,  well  within  the  5.7  W  uncertainty  calculated 
in  Appendix  B.  Figure  20  shows  a  sampling  of  power  settings 
and  the  impact  of  the  differences  in  specific  heat  values. 

The  possibility  of  contamination  of  the  coolant  was 
explored,  although  all  fluid  circulation  components  were 
flushed  with  the  coolant  prior  to  experiments.  It  was  noticed 
late  in  the  data  collection  process  that  the  tubing  in  the 
pump  head  had  split.  An  examination  of  the  tubing  and 
comparison  with  new  tubing  revealed  a  hardening  had  occurred 
during  the  course  of  usage,  indicating  the  likelihood  of  some 
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sort  of  reaction  between  the  Tygon  tubing  and  the 
polyalphaolefin.  A  conversation  with  the  coolant  manufacturer 
supported  this  as  a  possible  source  of  contamination.  The 
progressive  nature  of  such  contamination  would  expectedly 
cause  data  from  early  experiments  to  differ  from  data 
collected  later  at  the  same  settings,  however  this  was  not  the 
case. 

Due  to  the  variations  in  heater  resistances,  power  into 
each  heater  was  not  identical.  It  was  noted  that  the  power 
into  heaters  D  and  F  was  usually  lower  than  that  of  the  other 
heaters  by  approximately  0.1  W.  This  might  account  for  some 
minor  diminishing  of  the  temperature  rise  which  would  have 
occurred  if  all  heater  powers  were  uniform.  Wire  resistances 
between  multiplexers  and  lead  wires  were  measured  for  both 
heaters  and  precision  resistors  and  found  to  be  0.1  ohms  or 
less  in  all  cases.  Hewlett-Packard  manuals  indicated  that  the 
voltage  measurements  of  the  HP  44701A  voltmeter  such  as  was 
utilized  vary  only  by  +(0.008%  +  700E-6)  volts 
(Hewlett-Packard,  1987) .  Uncertainty  for  the  power  to  each 
heater  was  thus  very  small  (see  Appendix  B  calculations) . 

Given  the  results  and  conclusions  drawn  from  the 
subseguent  analysis,  it  is  suspected  that  the  measured 
temperature  differential  (outlet  minus  inlet  temperature)  was 
not  sufficiently  accurate  due  to  the  slight  fluctuations  in 
inlet  temperature,  particularly  at  higher  flow  rates.  At 
lower  flow  rates,  the  pump  produced  a  more  "pulsed”  flow  that 
may  have  resulted  in  degraded  accuracy  in  the  inlet  and  outlet 
temperature  measurements. 
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VI.  RECOMMENDATIONS 


In  order  to  minimize  uncertainty  in  further  research  with 
the  FTM,  three  recommendations  are  made.  First,  use  of  a 
larger  volume  refrigerated  bath  might  aid  in  stabilizing  inlet 
temperatures  by  allowing  greater  resident  time  for  coolant 
prior  to  its  entrance  to  the  module.  Additionally,  efforts 
should  be  made  to  place  thermocouples  in  the  inlet  and  outlet 
tubing  closer  to  the  actual  module  inlet  and  outlet.  The 
combination  of  these  two  changes  would  lend  greater 
credibility  to  the  resulting  energy  balance  by  reducing  the 
uncertainty  in  the  fluid  temperature  differential.  Finally, 
efforts  should  be  made  to  determine  the  uncertainty  of  the 
coolant's  thermophysical  properties  prior  to  further 
experimentation . 

A  closer  examination  of  the  surface  temperatures  in  the 
vicinity  of  the  inlet  channel  (as  shown  in  Figure  6)  through 
the  emplacement  of  additional  thermocouples  would  be  useful  in 
the  study  of  the  count erf low  heat  exchanger  effect  in  that 
region.  Thermocouple  placement  along  the  right  side  of  the 
module  in  the  areas  where  flow  changes  direction  and  exits  the 
FTM  would  aid  in  the  analysis  of  the  heat  transfer  mechanisms 
in  operation  there.  Analytical  modelling  of  the  FTM  could  be 
enhanced  by  data  obtained  through  these  changes. 

An  x-ray  of  the  module  would  provide  valuable  information 
regarding  the  internal  structure.  Of  particular  interest 
would  be  the  structure  in  the  areas  where  flow  direction 
changes  and  in  the  inlet  and  exit  channels. 

Finally,  the  use  of  thermochromic  liquid  crystals  would 
provide  a  complete  visual  image  of  the  surface  temperature 
patterns  at  various  power  and  flow  settings.  Pictures  of 
these  images  could  be  compared  with  the  temperature  data  for 
verification. 
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APPENDIX  A 


SAMPLE  CALCULATIONS 


Throughout  this  appendix,  the  sample  calculations  are 
based  upon  power  setting  10  at  a  flow  rate  of  1000  ml/min. 
Data  used  for  calculations  is  shown  in  Tables  17  and  27. 


1. 


Calculation  of  Coolant  Thermophysical  Properties 
From  Equation  4, 


T 

““■avg 


(  T(4)  T(3)  ) 

2 


19.37°C 


Density  was  calculated  using  Equation  5  as  follows: 
p  =  0.794  -  (  (8.5E-4)  *  (19.37  -  20)  )  =  0.7945  g/ml 


Specific  heat  was  determined  via  Equation  6  as: 

Cp  =  (0.52  +  (7.1E-4  *  (19.37  -  10)  )  *  4.184  =  2.204  J/(g°C) 


2.  Calculation  of  Heat  Dissipation 

Table  24  shows  the  data  used  in  Equation  7  to  determine 

ATg^g  =  2.61  as  obtained  by  calculation  within  the  data 

collection.  Utilizing  this  value,  heat  absorbed  by  the 

fluid  was  calculable  as: 

=  (0.7945)  ^  (2.204)  *  (2.61)  ^  ^ 


Insulation  losses  were  determined  using  Equation  9  as: 
=  (0.04)  *  ( 0  ._0,.^1_9  ).  (19.62  -  18.20)  =  0.07  W 

«ins  0.0191  '  ' 


Dividing  the  above  result  by  area  results  in: 

q"l„  =  -^7^  (19-62  -  18-20)  ■=  2.97  W/m" 
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Table  24.  Calculation  of  . 

civy 
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APPENDIX  B.  UNCERTAINTY  CALCULATIONS 


Uncertainty  calculations  were  evaluated  using  the  root 
mean  square  method  (Beckwith,  Marangoni,  and  Lienhard,  1993) . 
Table  24  provides  a  summary  of  the  uncertainty  in  variables 
used  for  calculations. 


VARIABLE 

UNITS 

UNCERTAINTY 

REFERENCE 

V 

ml/min 

maximum  of 

2.8% 

in  range  of 
testing 

Flowmeter  calibration 

(see  Table  26  and 
corresponding  curve 
fit  in  Figure  41) 

T 

°C 

0. 13°C 

Thermocouple 

calibration  (see 
Table  27) 

V 

volts 

0.008% 

+ 

700E-6 

Integrating  voltmeter 
manual  (Hewlett- 

Packard,  1987) 

R 

ohms 

0.0001  ohms 

Precision  resistor 

calibration 

A 

square 

meters 

10%  (0.00219 

square  meters) 

Estimated 

L 

meters 

0.003  m 

Ruler  resolution 

Table  25.  Uncertainties  in  variables. 

Uncertainty  in  AT^^g  =  (avg) -T^^(avg) )  was  calculated 

as: 


Uat,.,  =  ^  (Ut)"  +  (Ut)2  (10) 

Thus,  for  the  stated  uncertainty  in  T, 

U.T  =  (0.13)2  +  (0.13)2  =  0.18°C  (11) 

•“■•‘■avg 
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For  power  setting  10  (69.68  W  input)  and  a  flow  setting 
of  1000  ml/min 

p  =  0.7945  g/ml,  Cp  =  2.203  J/ (g°C)  , 

=  2.61°C,  ATin,  =  0.40°C, 

k  =  0.04  W/(m°C),  L  =  0.0191  m,  and 

A  =  0.0219  m^. 


The  uncertainty  in  was  calculated  as: 


Q,  =  )  ^  -  (  (  (V)  ) 


5.7  W 


Uncertainty  in  the  insulation  losses  is  calculated  as: 


Similarly  Ug//  is  determined  as  follows: 
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Uncertainty  in  the  power  calculations  was  determined  as: 


(15) 


For  heater  A  at  power  setting  10,  with 
Vjj  =  16.27  volts,  Vj,  =  0.35  volts,  and 

R  =  0.4876  ohms  (determined  as  discussed  in  Chapter  II), 
calculations  yield  Up^  =  0.07  W. 

Specific  heat  data  from  the  manufacturer  of  the 
polyalphaolef in  and  testing  results  performed  by  an 
independent  laboratory  are  plotted  in  Figure  42.  Testing  data 
was  extrapolated  to  show  values  in  the  temperature  range  of 
this  research. 
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%FLOW 

MIN/ (2L) 

UNCERTAINTY  (MIN) 

ML/MIN 

%UNCERTAINTY 

0 

0.00 

0.00 

0 

0.0 

10 

34.27 

0.17 

58 

0.0 

15 

18.00 

0.17 

111 

0.1 

20 

12.50 

0.17 

160 

0.2 

25 

8.00 

0.13 

250 

0.4 

30 

5.77 

0.10 

347 

0.6 

35 

4.63 

0.10 

432 

0.9 

40 

4.02 

0.10 

498 

1.2 

45 

3.53 

0.07 

566 

1.1 

50 

3.13 

0.07 

638 

1.4 

55 

2.87 

0.07 

698 

1.6 

60 

2.60 

0.07 

769 

2.0 

65 

2.37 

0.07 

845 

2.4 

70 

2.17 

0.05 

923 

2.1 

75 

2.05 

0.05 

976 

2.4 

80 

1.90 

0.05 

1053 

2.8 

85 

1.80 

0.05 

1111 

3.1 

90 

1.68 

0.03 

1188 

2.4 

100 

1.50 

0.03 

1333 

3.0 

Table  26.  Flowmeter  calibration  with  estimated  uncertainty. 


PLATINUM  RESISTANCE 

INLET  T(3) 

OUTLET  T(4) 

THERMOMETER  (°C) 

(°C) 

(°C) 

17.101 

16.983 

16.972 

17.841 

17.742 

17.721 

19.380 

19.271 

19.272 

25.149 

25.091 

25.078 

30.275 

30.197 

30.185 

30.294 

30.216 

30.180 

36.043 

36.044 

36.014 

41.980 

41.936 

41.925 

46.376 

46.354 

46.347 

Table  27.  Calibration  of  inlet  and  outlet  thermocouples. 
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Calibration  of  Floy  tleter 

6  Cc|n  183  y-(atbx+cx2+dx3f9x'l  +  fx5)  r  2=0.  999555092 
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Figure  41.  TableCurve  polynomial  fit  of  flowmeter 

calibration  data. 
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37.  5  62.  5  87.  5 

P0rcent  Flou 


Specific  Heat  [  cal/(g  *  deg  C)  ] 


Temperature  (degrees  Celsius) 


Figure  42.  Comparison  of  specific  heat  values. 
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